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ABSTRACT 

The cooling history of a quark star in the color superconductive phase is investigated. Here we 
specifically focus on the 2-flavour color (2SC) phase where novel process of photon genera- 
tion via glueball (GLB) decay have been already investigated (Ouyed & Sannino 2001). The 
picture we present here can in principle be generalized to quark stars entering a superconduc- 
tive phase where similar photon generation mechanisms are at play. As much as 10 45 — 10 47 
erg of energy is provided by the GLB decay in the 2SC phase. The generated photons slowly 
diffuse out of the quark star keeping it hot and radiating as a black-body (with possibly a 
Wien spectrum in gamma-rays) for millions of years. We discuss hot radio-quiet isolated neu- 
tron stars in our picture (such as RX J185635-3754 and RX J0720.4-3125) and argue that 
their nearly blackbody spectra (with a few broad features) and their remarkably tiny hydro- 
gen atmosphere are indications that these might be quark stars in the color superconductive 
phase where some sort of photon generation mechanism (reminiscent of the GLB decay) has 
taken place. Fits to observed data of cooling compact stars favor models with superconduc- 
tive gaps of A2SC ~ 15-35 MeV and densities P2SC = (2.5-3.0) x (pn being the nuclear 
matter saturation density) for quark matter in the 2SC phase. If correct, our model combined 
with more observations of isolated compact stars could provide vital information to studies of 
quark matter and its exotic phases. 
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1 INTRODUCTION 

Quark matter at very high density is expected to behave as a color 
superconductor (e.g. Rajagopal & Wilczek 2001). Of interest to 
the present work is the 2-flavor color superconductive phase, the 
2SC phase, where only the up and down quarks of two color (say, 
Ml, «2, di, and d2) are paired in a single condensate, while the ones 
of the third color (say, 113 and d$) and the strange quarks of all 
three colors are unpaired. Associated with superconductivity is the 
so-called gap energy A2SC inducing the quark-quark pairing and 
the critical temperature (T c ) above which thermal fluctuations will 
wash out the superconductive state. 

The cooling of quark stars in such a phase 1 has previously 
been investigated in the literature. It was found that a crust must be 
included for 2SC star cooling to be compatible with existing data 
(Blaschke et al. 2000). As for the cooling of hybrid stars (neutron 



1 Here we only consider the stars in 2SC phase. Cooling of quark stars 
in the Color-Flavor-Locked phase where the strange quark enters the dy- 
namics, has also been investigated in the literature and it was shown they 
cool down too rapidly in disagreement with observations (e.g., Schaab et al. 
1997; Blaschke et al. 2001). 



stars with a 2SC core), it is believed that small gaps (A2SC < 1 
MeV) tend to reproduce cooling curves which agree well with the 
observed data (we refer the interested reader to §6 in Weber 2005 
and references therein). Whether such small gaps can be justified 
and why the crust - believed to be tiny in quark stars - is so crucial 
to cooling is debatable. More recent studies of these hybrid stars 
find that unless all quarks are gaped the stars cool too fast in dis- 
agreement with observational data (the fast cooling is due to direct 
Urea process on unpaired quarks; we refer the reader to Blaschke 
et al. 2005 for more details). However the studies mentioned above 
did not include internal heating mechanisms as we do here (the 
GLB decay) which as we will show have definite consequences on 
the cooling history and spectral features of these stars. Furthermore, 
there are still pulsar candidates these models fail to explain leaving 
room for exploring other possibilities. 

A novel feature of the 2SC phase is the generation of GLB par- 
ticles (hadrons made of gluons) which as demonstrated in Ouyed & 
Sannino (2001) immediately decay into photons. These GLBs ap- 
pear at temperatures much below the critical temperature for the 
onset of the superconductive phase. We isolate three cooling steps: 
(i) in the early hot stages when the quarks are still unpaired the 
neutrino emission is due to the three quark direct Urea processes 
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u c + e~ —* d c + u e and d c —* u c + e~ + where the subscript c 
denotes the quark color. The star rapidly cools to T c thus entering 
the 2SC phase; (ii) in the 2SC phase the paired up and down quarks 
of two colors do not contribute to cooling thus reducing neutrino 
emission (see Page&Usov 2002 for more details). While reduced 
by a factor of 3, the neutrino emission remains the dominant cool- 
ing mechanism until the star has cooled to the critical temperature 
for GLB formation; (iii) at this point in time, the GLBs decay to 
photons providing an extra heat source and the cooling becomes 
driven by these photons slowly escaping from the star. We should 
note that this photon cooling stage might be particularly dominant 
for the further cooling history of the star for cases and situations 
where neutrino cooling is quenched (e.g., by secondary pairing of 
the «3 and quarks as discussed in Page&Usov 2002). This third 
phase, unique to our model, alters the cooling history of the star 
keeping hot and radiating as a blackbody for millions of years as 
we show in this work. 



1.1 Caveats on the 2SC phase and model generalities 

It has been argued that once the conditions of charge neutrality 
and /3-equilibrium are enforced it appears that the so-called "gap- 
less 2SC phase", or g2SC, might be the favoured ground state 
(Shovkovy & Huang 2003; Aguilera et al. 1995). Other studies in- 
stead favor the so-called "gapless CFL phase", or gCFL (Alford, 
Kouvaris, & Rajagopal 2004; and references therein). These alter- 
natives (g2SC and gCFL) however might be prone to instabilities 
begging for serious scrutiny before this debate can be settled (e.g. 
Shovkovy, Riister, & Rischke 2004; Huang & Shovkovy 2004; He 
et al. 2005). The differences in these studies are connected with the 
choice of the model parameters and is beyond the scope of this pa- 
per. For g-2SC matter the most likely resolution of the instability 
is a transition to a crystalline superconductor (e.g. Casalbuoni et al. 
2005) 2 . If such a transition occurs on much longer timescales than 
GLB formation and decay time (i.e. 10 -13 -10 -14 s; Ouyed & San- 
nino 2001) then our model should still be valid; a stable 2SC/g-2SC 
phase is not a necessary condition. As we have said, the picture we 
are presenting here in its generality should apply to quark stars in 
any phase where photon generation mechanisms are at play (e.g. 
Vogt, Rapp, & Ouyed 2004). 

The paper is presented as follows: In Sect. 2 we briefly present 
the GLB formation and decay in the 2SC phase followed in Sect. 
3 by a description of the model's assumption. Cooling calculation 
results are shown in Sect. 4 with a comparison to observed data. 
In particular, radio-quiet isolated neutron stars are discussed in the 
context of our model in Sect. 5. We conclude in Sect. 6. 



2 GLUEBALL FORMATION AND DECAY 

The GLB mass is given in Ouyed&Sannino (2001) as, 
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Figure 1. Glueball mass (solid) and deconfinement temperature (dotted) vs 
chemical potential for two different QCD energy scales, Aqcd = 150 
MeV (top panel) and Aqcd = 300 MeV (lower panel). 

where b = 22/3, e 2 /47r = 1/137 and d is a positive constant of 
order unity, and 



A ~ A exp 



/'• 



11 



(2) 



In the expression above, fi oc p 1 ^ 3 is the chemical potential (p is 
the quark matter density) and g(p) the effective energy dependent 
coupling constant define as p(p) 2 /4-7T ~ 1/ 1ii(.E/Aqcd) where 
Aqcd is the QCD energy scale. Here E is the energy scale at which 
the theory is being applied which in our case is E ~ p. 

It has been shown in Sannino et al. (2002) that the GLBs ap- 
pear (or enter the dynamics) at the critical temperature, 



Tglb = 
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(for T > Tglb, GLBs "melt") where v = 1/vAeis the gluon ve- 
locity with the dielectric constant (e) and the magnetic permeability 
(A) given as (Rischke et al. 2001) 
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Figure Q shows the mass, Mglb, and Tglb of the GLBs for 
different gaps and chemical potentials for the two cases of Aqcd = 
150 MeV and 300 MeV. For the quark matter densities considered 
here and a gap of tens of MeV the GLBs mass is of the order of 
a few MeV. We also note that in general Tglb — Mglb- In the 
following we define E-, as the photon energy from the GLB — > 77 
reaction. 



2 It remains to be shown that such a phase possesses a thermodynamic 
potential (free energy) which is below the unstable g-2SC phase. While this 
might be the case in a narrow density range in the QCD phase space, more 
studies involving more complicated crystalline structures seem to be needed 
to settle this issue (Casalbuoni 2003). 



3 SCENARIO DESCRIPTION AND ASSUMPTIONS 

We begin with a hot quark star in the quark-gluon plasma phase 
cooling and undergoing the transition into the 2SC phase. As men- 
tioned in the introduction, neutrino cooling will be dominant at first 
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and the star rapidly reaches the temperature T c ~ 0.57A2SC for 
the onset of the 2SC phase. The neutrino cooling is also dominant 
in the 2SC phase despite a decrease in neutrino emission induced 
by pairing of quarks (Carter&Reddy 2000; Page&Usov 2002) until 
ultimately the star cools to Tglb at which point in time the GLBs 
start forming. We assume for simplicity that the entire star cools 
to Tglb and any temperature gradients in the early stages will be 
smoothed out; which is not unreasonable given the extremely large 
conductivity of the 2SC phase (e.g. Heiselberg&Pethick 1993). All 
of the GLBs immediately decay into photons and beyond this stage 
cooling is dominated by the photons which slowly diffuse out of 
the star. 



3.1 Total photon energy and plasma cut-off 

It was shown that 3/8 of the total number of gluons will form 
GLBs (Ouyed & Sannino, 2001). This implies that 3/8 of the nu- 
cleon binding energy, 938 MeV, is transformed to GLBs which 
amounts to roughly 350 MeV in photon energy per nucleon. For 
a 10 km radius quark star and matter densities on the order of 2/9n 



(p ~ 350 MeV) this translates to a total energy E 1 



10° 



ergs released as photons; pN = 2.71 x 10 14 g cm" 3 is the nu- 
clear matter saturation density. However, let us recall that the emis- 
sivity of photons at energies below that of the plasma frequency 
(E p = Tiuip ~20-25 MeV for the density regime representative 
of 2SC phase) is strongly suppressed (Alcock, Farhi, & Olinto 
1986; Chmaj, Haensel, & Slominsli 1991; Usov 2001). As such, 
only the tail (E-y > 20 — 25 MeV) of the distribution will ex- 
ist inside the 2SC star. The true available total energy in photons 

isthen£ 7 , t ~ [/~/ 7 d£ 7 // °°/ 7 d.E 7 l x £ T , to t. ^ (HT 6 - 

10" 8 ) x £ 7 ,tot. (10 45 -10 47 ) ergs; here / 7 is the Planck distri- 
bution; we note that different distributions lead to same values for 
.E 7 ,t but as we state in §4.1 these photons will eventually acquire a 
Wien distribution as they leak out of the star. 



4 COOLING AND LIGHTCURVES 

For the range in temperature, T, considered here the vast major- 
ity of quarks in the star are in their ground states kT/Ep « 1 
where Ep is the Fermi energy. The density of the free quarks at 
temperature T that will interact with the photons is 



N a = - 
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2h 2 



TT 2 h 
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where e is the electronic charge and m q is the quark mass. 

The problem we are solving here is that of a uniform sphere 
of constant density and constant temperature gas of fermions and a 
homogeneous source of photons with E-, > fcT q where T q is the 
the temperature of the quarks in this case. The number of photons 
escaping as a function of time we calculated using the numerical 
method described by Sunyaev&Titarchuk (1980, hereafter ST80; 
see also Shapiro et al. 1976 and Miyamoto 1978). The correspond- 
ing diffusion equation in spherical coordinates is 
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where J is the average intensity of emission and is a dimensionless 
quantity (following ST80 choice of units), u is the dimensionless 
time (u = caN q t; c is the speed of light, t is the time in seconds, 
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Figure 2. The total energy that has escaped from the surface of a 10 km ra- 
dius quark star is plotted for increasing time. The amount of escaped energy 
is shown as a percentage of the initial energy contained within the quark 
star at t = 0; the curves from left (Tqlb = 1 MeV) to right (T GLB = 10 
MeV) are in 1 MeV step. 



is the quark wavelength ), and r is the dimensionless optical dis- 
tance from the center of the star (r = aN q r where r is the distance 
from the center of the star). To find the escape time distribution 
function P(u) of a single photon we have to normalize the func- 
tion obtained J (to, u) (see Appendix A in ST80 for details) so that 
f °° P(u)du = 1, and 
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exp 
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where to = <rN q R q is estimated at surface of the star of radius 
R q . Equation gives the normalized probability of escape for 
photons from the surface of the star as a function of time. We then 
adjust the magnitude of the distribution so that its area is equal to 
the total number of photons, with i? 7 > E p , from GLB decay. The 
total number of photons escaping versus time was then converted 
to luminosity and to surface temperature, T a , using the black-body 
formulation (see below). 

Figure |2|shows the cumulative percentage of the total energy 
generated through GLB decay that has been radiated from the sur- 
face of 10 km radius star as a function of time. It takes at least 
5 x 10 5 years for the energy to escape the dense quark star and 



7r X q is introduced to make use of dimensionless 



and a = 7rA q where A q 



(2Trhc/^/(k B T) 2 + m 2 c 4 ) 



3 The cross-section a 
description as explained in ST80 following Chapline&Stevens (1973). We 
should note however that a complete mechanical prescription (using rela- 
tivistic Compton cross-section in the appropriate regime where hv > kTq) 
was used in deriving equation^] More specifically the treatment of cross- 
section we adopted is that of Miyamoto (1978) which is consistent with the 
use of the Compton Focker-Planck equation. 
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so we conclude that GLB decay in quark stars does not constitute 
a suitable inner engine for gamma ray bursts (GRBs) as originally 
claimed in Ouyed&Sannino (2002); recent work by Ouyed et al. 
(2005) suggests instead that quark stars in the CFL phase consti- 
tute better candidates for GRB inner engines. However, Fig. 2 and 
Fig. 3 show that if GLB decay does occur in a quark star the object 
will be extremely bright and should be readily observable. 

4.1 The escaping photons: A Wien spectrum 

In its simplest approach the quark-photon interaction can be de- 
scribed using elementary random-walk theory which tells us that 
the average number of quark scatterings will be approximately 
equal to t s 2 , with r s = R/\ v and A„ ~ 1 fm being the Comp- 
ton mean-free path of the photons. Furthermore, the average frac- 
tional energy transfer in a quark-photon scattering will be of the 
order of (2kT q /m q c 2 ). For r s >> 1, as it is in our case, one can 
show that the initial photon spectrum is distorted into a distribution 
having approximately the shape of a Wien spectrum (e.g., Chapline 
& Stevens 1973). In general, it is found that an equilibrium so- 
lution to the Compton Fokker-Planck equation is a Bose-Einstein 
distribution which, for hv > fcT q and large r, approaches a Wien 
spectrum (Miyamoto 1976 and references therein). That is, by the 
time the photon escape the star their distribution will be given by 
F v oc i; 3 exp(-/u//fcT). 

4.2 Thermalization: A black-body spectrum 

A property of quark stars is the existence of an "electrosphere" ex- 
tending a few hundred to a thousand Fermi above the surface of 
the star (see Usov, Harko & Cheng for a recent study) 4 . Based on 
a collisional treatment of photons on single electrons and given the 
electron density in the electrosphere, n e ~ 10~ 5 — 10~ 4 fm -3 , one 
can show that photons with energies i? 7 > 6.5MeV will not be af- 
fected by the layer (Cheng & Harko 2003; Jaikumar et al. 2004). In 
other words, the escaping photons do not scatter often in the elec- 
tropshere to become thermalized. 

However, the outflowing high energy (E~, >> m c c 2 ) pho- 
tons enter a region rich in (e + e~) pairs and low energy (<< E y ) 
photons (Aksenov et al. 2004). The corresponding 'compactness 
parameter' I = £ 7 <7t/47t i?c 2 >> 1 (e.g. Frank, King, & Raine 
1992) implies ideal conditions for 77 — > e + e~ processes to take 
place quickly leading to thermal equilibrium (e.g. Paczyhski 1990). 
We thus expect the Wien distribution of the escaping photons to 
evolve into a black-body although a residual Wien tail in gamma- 
rays (Js 20-25 MeV) could survive as a signature of the original 
distribution. 

4.3 Crust effects 

The electrosphere allow for radial electric fields whose magni- 
tude can support matter against the star's gravity. The mass of the 
crust which can exist suspended by this electric field cannot exceed 
10" 5 M Q (Alcock et al. 1986; Horvath et al. 1991). Its effect can be 
quantified by relating the "internal" temperature T s or the temper- 
ature at the bottom of the crust to the effective temperature at the 

4 A charge neutralizing surface electron layer can be realized in the 2SC, 
the gapless CFL (gCFL) or crystalline color superconductive phases (Al- 
ford&Rajagopal 2002; Shovkovy&Huang 2003), but not in the pure CFL 
phase (see however Usov 2004). 



surface of the crust via the so-called Tsuruta law (Tsuruta 1979). 
The formula for thick crusts is (see also Shapiro&Teukolsky 1983, 
p330) T c = (10T S ) 2/3 . We note however that for quark stars the 
crust is much thinner than in the NSs. Therefore, we will investigate 
the two extreme cases of negligible or no crust (T c ~ T s ; e.g. Piz- 
zochero 1991) and thick or maximum crust (M max . ~ 10 M©) 
using the Tsuruta formula. In both cases, the effective temperature 
as seen by a distant observer is T c °° = T c \J\ — Rsch./R where 
^Sch. is the star's Schwarschild radius. 

4.4 Comparison to observations 

The related cooling curves are shown in Figure [5] for different 
initial temperatures To and for a star of radius 10 km. Models 
with maximum crusts do not appear to be successful when com- 
pared to cooling of neutron stars if these are quark stars. The cor- 
responding cooling curves start overlapping with observed data 
for Tglb > 0.57A2SC which is outside the 2SC phase where 
GLBs cannot form. Such temperatures will not be sustained since 
we expect neutrinos will drive Tglb to lower values early in the 
star's cooling history. When compared to observed data our cool- 
ing curves agree best with models of stars with thin or no crusts 
and for temperatures 0.8 < Tglb < 2 MeV. The regime favored 
by the best fits correspond to densities p2SC = (2.5-3.0) x pn (or 
^ = 350-450 MeV) and energy gaps A 2S c = 15-30 MeV (§2). 

In general, in our model best fits to observations favor quark 
stars with thin or no crusts. This is also consistent with the fact that 
quark stars are expected to be born bare given the extreme tem- 
peratures involved during their formation. In the so-called 'Quark- 
Nova' scenario (Ouyed et al. 2002; Keranen et al. 2005) for exam- 
ple the crust could be regenerated from fall-back material but it is 
expected to be tiny. This brings to the special family of hot radio- 
quiet isolated neutron stars (INSs) within the picture we presented 
so far. Specifics below. 



5 RADIO-QUIET ISOLATED NEUTRON STARS 
5.1 Observations 

There are about 7 radio-quiet isolated neutron stars (INSs) that 
show features that might be of interest to us. There are important 
clues from observations that are worth mentioning: (i) While re- 
cent XMM-Newton observations of 4 radio-quiet INSs show clear 
features (in particular RX J0720.4-3125 show a clear broad feature 
at around 310 eV; Haberl et al. 2004), RX J 1856.5-3754 remains 
featureless; (ii) it is believed that the observed broad features are 
absorption due to the proton cyclotron line. These hydrogen lines it 
was suggested could arise from a pure hydrogen atmosphere (van 
Kerkwijk et al. 2004); (iii) what is further remarkable about these 
objects is the optical excess observed in their spectra as compared 
to the X-ray component (e.g. kaplan et al. 2003). Recent studies 
of this phenomenon show that a thin (1 cm thin and not exceeding 
10 14 g) layer of hydrogen atmosphere can well represent both the 
optical and X-ray data without invoking additional thermal com- 
ponent (Motch et al. 2003); (iv) finally, there are indications that 
accretion might not be important since the rates must be many or- 
ders of magnitudes smaller than the Bondi-Hoyle rates of ~ 10 s 
g/s (Bondi & Hoyle 1944; Bondi 1952). Motch et al. (2003) esti- 
mate from fits to optical emission that the total mass of hydrogen 
is no more than 2 x 10 12 g which translates into accretion rates not 
exceeding 1 g/s. A second clue that accretion might not be at play 
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Figure 3. The time evolution of the temperature, T c °° (as seen by a distant observer) for a star of a 10 km radius. Shown are cooling curves for different initial 
temperatures Tqlb- Solid lines are for a star with no crust while the dotted ones is for the maximum crust case. The cooling curves begin approximately 2 
weeks after t = 0. Experimental data taken by ROSAT and ASCA (Table 7 in Weber 2005) are also shown. 



in these sources is the H a nebula surrounding RX J1856.5 — 3754 
which can be attributed to a relativistic wind (van Kerkwijk & 
Kulkarni 2001). 



5.2 Our scenario 

The newly born bare quark star could regain its crust by accre- 
tion from the ISM or from fall-back material during its formation 
(Keranen et al. 2005). In both of theses possibilities a tiny hydrogen 
envelope can only be explained if accretion rates are less than 1 g/s. 



Such rates are puzzling and difficult to account for and we would 
rather assume that accretion is irrelevant (see point (iii) in §5.1) and 
that the star regained a crust from fall-back material as in the case 
of the 'Quark-Nova' scenario (Ouyed et al. 2002; Keraanen et al. 
2005) or any similar formation scenarios of quark stars. We should 
however mention the scenario discussed in Usov (1997) where it 
is argued that in the case of more or less spherical (Bondi) accre- 
tion onto a nearly bare quark star only a small part of falling gas is 
collected at the stellar surface. The bulk of this gas passes through 
the quark surface leaving a rather low-mass steady hydrogen atmo- 
sphere (~ 10 12 g). Below we discuss another possibility whereas 
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the original crust is slowly depleted in time that might possibly ex- 
plain the clustering in time and temperature of INSs. 

An important intrinsic property of quark stars is the depen- 
dency of the size of the crust on the properties of the electropshere 
(Alcock et al. 1986; Usov et al. 2005) which in turn are very sen- 
sitive to the global charge neutrality of the star and its temperature. 
While in the case of pure CFL it has been argued in the literature 
that in such a state the star is globally neutral by construction (Ra- 
jagopal & Wilczek 2001 ; see however Usov 2005) we ask ourselves 
the question of what happens in the case of 2SC, g2SC and gCFL? 
As these stars cool one could imagine that they slowly evolve to a 
more globally neutral configuration. Presumably the electrosphere 
and the corresponding electic field are reduced in time which would 
correspond to the crust shrinking to eventually reach infinitesimal 
size when most of the electron in the electropshere have been reab- 
sorbed by the quark matter. What is interesting in this scenario we 
suggest - which remains to be confirmed - is the fact that the bottom 
of the crust (dominated by the settled heavy elements) sinks closer 
to the quark matter surface. Over time the crust material is slowly 
depleted/deconfined from the bottom leaving the lighter upper at- 
mosphere material. 

In our scenario then the smaller the crust gets as the star cools 
and reduces its electrosphere the richer in lighter elements (and ulti- 
mately in hydrogen) it becomes. The end result we suggest is a tiny 
envelope reminiscent of the thin hydrogen envelope inferred from 
observations. Hence we can define a universal critical temperature 
Th below which the envelope becomes transparent leading to the 
optical excess. The rate of depletion of the envelope closely fol- 
lows the cooling rate/curve. As such, the size of the envelope will 
drop quicker as the stars enter the faster cooling beyond the 10 5 - 
10 6 years epoch. This one-to-one relationship could in principle 
account for the observed clustering in time and temperature of the 
7 radio-quiet INSs as well as their their associated optical excess. 
The critical temperature Th among others features inherent to our 
model and the observational consequences will be explored in more 
details in an upcoming paper. For now we are tempted to specu- 
late that radio-quiet INSs are quark stars in the superconductive 
phase cooling through the slow release of photons from GLB de- 
cay or similar processes. For completeness, shown in Figure|4]are 
RX J185635-3754 and RX J0720.4-3125 (Haberl et al. 1997; Wal- 
ter&Lattimer 2002; Kaplan et al. 2003a; Walter et al. 2004) which 
line up closely with the Tglb = 1 MeV curve for a 10 km quark 
star. 




Time (years) 

Figure 4. The time evolution of the temperature, T c °° (as seen by a distant 
observer) for Tqlb = 1 MeV, and for different star radius R. The cooling 
curves begin approximately 2 weeks after t = 0. The observed data are RX 
J 185635-3754 (lower square) and RX J0720.4-3125 (upper square). 

erg of photon energy can be provided by the GLB decay radiated 
away by the star over millions of years. The persistent blackbody 
emission offers a different picture from other models of cooling of 
quark stars where the photon spectrum is predicted to vary in time. 
Our fit to cooling data of observed INSs suggests that these are 
quark stars with extremely thin or no crust where GLB decay has 
taken effect. The fits favor gap energy A2SC = 15-35 MeV with 
corresponding densities P2SC = (2.5-3.0) x this might be of 
relevance to Quantum-Chromodynamics. Let us state again that in 
principle our picture applies to quark stars where photon generation 
mechanisms are at play. If correct, our model combined with more 
observations of cooling of isolated compact stars, in particular hot 
radio-quiet ones, could be used to constrain other fundamental pa- 
rameters of quark matter and its exotic phases. 



6 CONCLUSION 

In this work we investigated the effect of photons from GLB de- 
cay on the cooling history of quarks stars in the 2SC phase. While 
the early cooling history is dominated by neutrino emission, we 
find that further cooling is driven by photons, from the GLB decay, 
slowly diffusing out of the star. This cooling mechanism could be 
potentially important in situations or cases where the neutrino emis- 
sion is heavily suppressed by pairing of quarks. Up to 10 45 — 10 47 
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5 The possibility that RX J1856. 5-3754 is a bare quark star has been put 
forward in the literature (Xu 2003; Turolla et al. 2004). The novelty in our 
model resides in the GLB decay as an extra energy source. It has also been 
suggested that this object might be a naked neutron star where cohesive 
effects induced by moderately strong magnetic field leads to condensation 
in the atmosphere and the surface external layers (Perez- Azorin et al. 2005; 
van Adelsberg et al. 2004). 
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